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INTRODUCTION 


A  computer  simulation  was  developed  for  a  complete  artillery  safing  and 
arming  (S&A)  mechanism  containing  a  straight-sided  verge  runaway  escapement. 
Figures  1  and  2  show  the  essentials  of  two  different  configurations  of  such  a 
mechanism  which,  in  addition  to  the  escapement,  contains  a  spin  driven  rotor  and 
a  two-pass  involute  gear  train. 

While  the  mathematical  model  of  the  verge  runaway  escapement  represents  a 
new  study,  the  description  of  the  rotor  and  the  gear  train,  together  with  the 
essence  of  the  computer  program,  is  based  to  a  considerable  extent  on  reference 
1,  which  deals  with  the  computer  simulation  of  an  S&A  mechanism  containing  a  pin 
pallet  runaway  escapement. 

As  in  reference  1,  the  three  motion  regimes  of  the  runaway  escapement  (cou¬ 
pled  motion,  free  motion,  and  Impact)  are  considered.  Again,  the  pallet  of  the 
escapement  may  have  an  arbitrarily-located  center  of  mass. 

Complete  derivations  are  given  for  the  mathematical  model,  Including  the 
configuration  variations.  A  detailed  description  and  operational  instructions 
are  given  for  the  associated  computer  program  whose  output  furnishes  instanta¬ 
neous  kinematic  conditions  of  the  escapement  components  and  the  number  of  turns 
to  arm,  as  well  as  all  non-impact  contact  forces.  This  program  was  extensively 
tested  for  the  configuration  no.  2  data  (fig.  2)  of  the  M739  fuze  S&A  mechanism. 
This  mechanism  simulation  required  approximately  26  turns  to  arm  when  a  system 
coefficient  of  friction  of  0.1  was  used. 


DESCRIPTION  OF  COMPUTER  PROGRAM 


The  following  gives  the  essential  steps  of  program  SANDA2V  which  is  listed 
in  appendix  A.  The  mathematical  model  of  the  verge  runaway  escapement  naturally 
differs  from  that  of  the  pin  pallet  escapement  as  shown  in  reference  1,  espe¬ 
cially  since  the  differential  equations  for  entrance-  and  exit-coupled  motions  in 
the  verge  escapement  are  not  the  same.  However,  the  programming  schemes  are 
essentially  Identical.  Thus,  with  the  exception  of  controls  pertaining  to  the 
aforementioned  entrance-  and  exit-coupled  motions,  the  program  runs  parallel  to 
the  one  associated  with  the  pin  pallet  escapment. 

The  following  outlines  the  essential  features  of  the  program  without  giving 
specific  control  details. 

The  main  program  starts  with  reading  in  and  writing  of  all  relevant  physical 
data.  This  is  followed  by  the  computation  of  gear  ratios,  fuze  body  angles, 
centrifugal  forces,  gear  train  constants,  and  earliest  and  latest  possible  values 
of  the  gear  angles  by  way  of  subroutine  ALFA,  as  well  as  the  initialization  of 
the  gear  angles.  The  simulation  begins  with  the  entrance-coupled  motion  at  a 
starting  angle  PHID,  which  represents  that  angle  $  of  the  escape  wheel  which  is 
associated  with  the  approximate  center  of  the  entrance  working  surface  of  the 
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pallet.  This  angle  then  corresponds  to  a  cumulative  escape  wheel  angle  PHITOT  of 

0". 


Coupled  Motion  (Location  1) 


Depending  on  whether  entrance-  or  exit-coupled  motion  takes  place,  differen¬ 
tial  equation  B-100  or  equation  B-H9  of  appendix  B  must  be  solved.  To  this  end, 
the  main  program  calls  on  an  available  fourth-order  Runge-Kutta  routine.1 2  The 
main  purpose  of  the  subroutine  FCT  is  to  present  the  second-order  differential 
equation  in  terms  of  two  first-order  opes  to  RKGS.  PHI(l)  and  PHI(2)  represent 
the  angle  tp  and  the  angular  velocity  <j>  ,  respectively.  The  computation  of  all 
parameters  of  the  differential  equation  takes  place  in  subroutine  FCT  as  well  as 
in  subroutines  IN2  and  KINEM. 

The  subroutine  KINEM  computes  current  values  of  g,  \p  and  i|>  ,  as  shown  in 
appendix  C,  as  well  as  of  the  moment  arms  Aj,  B|,  Cj,  and  Dj  as  shown  in  appendix 
D. 

Subroutine  IN2  computes  various  gear  mesh  parameters  and  instantaneous  mesh 
contact  angles,  as  well  as  the  signum  functions  s^,  s*,  s^,  and  s^.  In  addition, 
the  parameters  Aj  to  A^,  listed  in  appendix  B,  are  obtained.  It  is  to  be  noted 
that  the  parameters  Ap  A5,  Ag ,  Au>  an(1  depend  on  whether  entrance-  or 

exit-coupled  motion  is  involved.  The  parameters  A* 5  to  A^q  are  computed  in 
subroutine  FCT.  The  reason  for  this  division  is  found  in  the  limited  number  of 
arguments  permitted  in  any  one  subroutine.  Finally,  the  gear-indexing  operation, 
as  outlined  in  reference  2,  is  performed  with  the  help  of  the  angle  <f>. 

The  following  discusses  the  manner  in  which  the  Instantaneous  rotor  angle 
<t>  +  N  A  of  the  coupled  motion  differential  equations  must  be  expressed  in 

suDroutIneT  FCT.  Recall  that  is  the  initial  rotor  angle,  <f>T  is  the  total 

angle  of  rotation  of  the  escape  wheel,  and  Njj  stands  for  the  gear  ratio  between 

the  rotor  and  the  escape  wheel.  Since  the  angle  <p  with  the  Runge-Kutta  variable 

PHI(l)  varies  between  approximately  134°  and  144°  during  entrance-  and  between 
209°  and  216°  during  exit-coupled  motion,  the  total  escape  wheel  angle  <f>  can 
only  be  obtained  by  continuously  adding  the  increments  due  to  each  cycle  of 

Runge-Kutta  computations.  Thus, 

<f>T  "  ♦xot  +  ^ 


1  RKGS  Routine,  IBM  System/360  Scientific  Subroutine  Package  (360A-CM-0X3) ,  Ver¬ 
sion  III. 

2 

The  program  uses  the  symbols  AA1 ,  etc.  throughout.  This  should  not  be  confused 
with  the  symbols  AAj  to  AA.g,  which  are  used  in  the  combined  exit-coupled 
motion  differential  equation  (8-119). 
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where 


<J>  -  total  escape  wheel  angle  up  to  a  certain  Runge-Kutta  cycle. 

(This  is  represented  by  PHITOT  in  the  program.) 

-  increment  of  escape  wheel  during  this  Runge-Kutta  cycle. 

The  increment  A<j>  is  calculated  as  the  difference  between  the  latest  value  of 
PHI(l)  and  its  previous  one  which  has  been  stored  as  PHIPR.  In  this  manner, 
equation  1  becomes 


4>t  -  PHITOT  +  PHI(i)  -  PHIPR 


(2) 


Subroutine  FCT  also  decides  on  the  values  of  Ip»  and  I,R  as  required  by 
equations  B73  and  B74  as  well  as  equations  B150  and  Bill  (app  B) .  The  associated 
conditional  statements  assign  the  larger  values  for  these  combined  moments  of 
inertia  whenever  the  product  of  the  angular  velocity  and  the  angular  acceleration 
is  positive;  l.e.,  both  quantities  have  the  same  sign. 

The  associated  suboutlne  OUTP  is  responsible  for  printing  out  the  results 
•  "  • 

<p,  and  <J>,  together  with  the  current  values  of  time,  g,  i|>,  i|>,  and  PHITOT. 
Further,  all  coupled  motion  contact  forces  are  calculated  according  to  equations 
B183  and  B186,  and  the  maximum  values  of  these  forces  during  one-rotor  cycle  are 
determined. 


Free  Motion  (Location  5) 


The  differential  equations  of  free  motion,  as  given  by  equation  B187  for  the 
pallet  and  equation  B193  for  the  combined  escape  wheel,  gear  train,  and  rotor  are 
again  solved  by  the  Runge-Kutta  routine.  To  obtain  the  magnitudes  of  the  vari¬ 
ables  <|>  and  ijj,  as  well  as  their  derivatives  at  identical  times,  the  two  indepen¬ 
dent  second-order  differential  equations  are  transformed  into  four  simultaneous 
first-order  ones.  (While  only  the  two  first-order  equations  associated  with  each 
of  the  two  variables  are  actually  coupled,  the  routine  treats  all  four  as  if  they 
were  coupled  and  thus  produces  solutions  for  identical  time  Increments.)  These 
four  expressions,  which  are  presented  in  subroutine  FCTF,  are  of  the  following 
forms 


DX(1)  -  X(2)  s  (-  <fr) 

(3) 

DX(3)  -  X(4)  :  (-  |) 

(4) 

DX(2)  -  I-A56(X(2))2  +  A57X(2)  +  A5g 

55 

+A59sin«>lc  +  N31  (*t  +  X(l)  -  PHIPR))]  :  («<(.) 
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DX(4) 


7T~~  l~A21  "  Aj2X(4)  - 

51 

+  Aj^sin(Yp  -  ♦  -  *C)J  !  (-'10 


(6) 


The  subroutine  FCTF  also  computes  the  parameters  A^  to  A^g  and  calls  on  sub¬ 
routine  IN2  for  the  computation  of  all  gear-related  parameters.3 

The  associated  subroutine  OUTPF  computes  the  free  motion  contact  forces 
according  to  equations  B199  and  B200  and  finds  their  maxima  after  determining  the 
parameters  A^  to  Aeg.  The  computation  of  these  latter  variables  also  requires  a 
call  on  subroutine  TN2.  In  addition  to  the  above , #  a  continuous  count  of  PHITOT 
is  provided  in  OUTPF.  This  angle,  time,  <f>,  <j>,  \p,  ip,  and  the  contact  forces  are 
printed  out.  This  subroutine  also  makes  the  decision  whether  or  not  to  remain  In 
free  motion.  The  sensing  variables  f  and  g'  *  GP  (equations  E4  and  E5  of  appen¬ 
dix  E)  are  used  for  this  purpose. 


Impact  (Location  15) 


• 

The  subroutine  IMPACT  uses  the  current  values  of  the  angular  velocities  ^ 
and  4*.  and  computes  the  post  Impact  angular  velocities  p.  and  applying  equa¬ 
tions  1B211  and  B212  of  appendix  B.  (Note  that  the  moment  or  inertia  of  the 
escape  wheel  is  now  expressed  according  to  equation  B213  (app  B),  which  refers 
the  rotor  as  well  as  the  gear  train  Inertias  to  the  escape  wheel.) 


Reversal  of  Gear  Train  Motion  Due  to  Impact 


If  the  impact  torque  on  the  escape  wheel  is  sufficiently  large,  the  motion 
of  the  gear  train  may  be  temporarily  reversed;  l.e.,  the  escape  wheel  angular 
velocity  £  may  become  negative.  This  would  cause  the  friction  forces  between  the 
gear  teeth  and  at  the  various  gear  pivots  to  be  reversed  in  direction.  (The 
normal  forces  between  the  gear  teeth  remain  unaffected,  and  the  normal  bearing 
forces  are  obtained  in  the  usual  manner.)  This  change  in  the  direction  of  the 
friction  forces  is  expressed  for  both  coupled  and  free  motion  by  letting  the 
coefficient  of  friction  m  of  all  gear  train  components  become  negative  (discussed 
in  app  E  of  ref  1).  This  # is  accomplished  in  subroutine  IN2  by  the  following  use 
of  the  signum  function  $/|$|s 


MU  -  ABS(MU)  *  W | l | 


(7) 


3  Whenever  Ae.  ■  IpR  ■  0,  the  simulation  stops  because  of  the  division  by  zero. 
Should  this  occur,  FCTF  printB  "IPR  EQUALS  ZERO  -  SIMULATION  TERMINATED." 
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(The  coefficient  of  friction  associated  with  the  escapement  interface  and  the 
pallet  pivot  is  called  i>  and  is  read  into  the  program  as  MU1.)  Any  motion 
reversal  at  these  surfaces  is  accounted  for  by  the  signum  functions  s^  and  85, 
respectively. 

Termination  of  Computations 

Computations  are  terminated  whenever  the  geared  motion  of  the  rotor  ends. 
This  corresponds  to  «  PH1CUTD.  The  duration  of  the  subsequent  unretarded 
motion  of  the  rotor  is  assumed  to  be  negligible. 


COMPUTER  SIMULATION  OP  EXAMPLE  MECHANISM 


This  section  discusses  a  computer  simulation  of  a  modified  S&A  device  of  the 
M739  fuze.  The  mechanism  has  configuration  no.  2,  as  shown  in  figure  2,  and 
contains  a  newly-designed  involute  gear  train.  While  this  gear  train  has  the 
same  gear  ratio  and  individual  center  distances  sb  the  original  design,  each  of 
the  meshes  now  has  unity  contact  ratio.  The  simulation  of  this  mechanism  was 
accomplished  with  the  help  of  computer  program  SANDA2V  which,  with  a  sample  out¬ 
put,  is  listed  in  appendix  A.  It  was  run  for  30,000  rpm  to  obtain  maximum  con¬ 
tact  forces. 

The  following  shows  the  input  requirements  of  the  program,  explains  the 
various  output  data,  and  discusses  the  manner  in  which  the  "number-of-turns-to- 
arm"  l.s  obtained  for  a  given  spin  velocity. 


Input  Data  (Appendix  A) 


The  first  portion  of  the  output  repeats  all  input  data,  which  represent  the 
mechanism  parameters  of  the  M739  fuze.  These  are  listed  both  as  computer  vari¬ 
ables  and  as  symbols,  according  to  appendixes  C  and  F  of  this  report: 


*  Both  meshes  were  designed  with  the  help  of  computer  programs  INVOL11  and 
GEARPAKAM2,  originally  shown  in  Progress  Report  No.  11  of  the  "Development  of 
Automated  Design  Optimization  Technique  for  Safety  and  Arming  Devices”  (Con¬ 
tract  No.  DAAK10-79-C-Q251,  January  15,  1981).  Copies  of  this  report  may  be 
obtained  from  either  F.  R.  Tapper,  ARRADCOM  or  G»  G.  Lowen,  The  City  College  of 
New  York. 
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Escapement  Paraoatara 


A  ■  a  ■  0.226  in.  (5.740  mm)  ■  distance  between  Divots  0_  and  0. 

(fig.  2)  P  8 

B  ■  b  ■  0.168  in.  (4.267  mm)  ■  escape  wheel  radius 

C  ■  c  ■  0.13138  in.  (3.337  mm)  ■  pallet  radius  as  defined  by  figure 

P-1  of  appendix  F 

ALPHEN  *  aan  -  43.5352"  ■  entrance  working  surface  angle 

ALPHEX  -  aex  -  29.2981°  -  exit  working  surface  angle 

NT  -  4  *  number  of  escape  wheel  teeth  spanned  by  verge 

CONFIG  -  2  ■  Configuration  no.  2  (See  section  on  fuse  body  configuration 

no.  2  in  appendix  B.) 

EREST  *  er  ■  0  ■  Coefficient  of  restitution 

LAMBDA  ■  A  ■  92.93°  -  angle  between  entrance  and  exit  pallet  radii 

(figure  F-l  in  appendix  F) 

N  ■  22  ■  Number  of  escape  wheel  teeth 

For  details  of  the  above  nomenclature,  see  appendixes  C,  E,  and  F. 


Mass  Parameters  of  Components 

Ml  ■  mj  *  0.3165  x  10  ^  lb-sec^/in.  (5.552  x  10~^  kg)  -  mass  of  rotor 

M2  ■  m2  «  0.3275  x  10  ^  lb-sec^/in.  (5.745  x  10-^  kg)  ■  mass  of  gear  and 

pinion  no.  2 

M3  ■  m^  ■  0.2631  x  10  ^  lb-sec^/in.  (4.615  x  10*"^  kg)  ■  mass  of  escape 

wheel  and  pinion  no.  3 

MP  **  nip  ■  0.1640  x  10  5  lb-sec^/in.  (2.877  x  10  ^  kg)  *  mass  of  pallet 

11  ■  Ij  ■  0.1967  x  10~5  in. -lb-sec^  (2.226  x  10~7  kg-ra^)  -  moment  of 

inertia  of  rotor 

12  ■  I2  ■  0.4026  x  10"7  in. -lb-sec^  (4.556  x  10”^  kg-m^)  ■  moment  of 

inertia  of  gear  and  pinion  no.  2 

13  ■  !•}  *  0.2125  x  10-7  in. -lb-sec^  (2.405  x  10"^  kg-m^)  ■  moment  of 

inertia  of  escape  wheel  and  pinion 
no.  3 
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IP  *  I  -  0.1950  x  10" 7  in.-lb-sac^  (2.207  x  10”®  kg-m^)  *  moment  of 
p  Inertia  of  pallet 


General  Parameters 


RC1  -  rcl  -  0.0576  in.  (1.463  mm)  ■  distance  from  pivot  of  rotor  to  its 

center  of  mass 

RCP  »  rcp  ■  0.  ■  pallet  eccentricity 

RHOP  -  p  -  0.0227  in.  (0.577  mm)  -  pallet  pivot  radius 

RPM  -  30,000  -  spin  rate 

PHI ICO  -  -  -*120.134°  -  rotor  angle  in  starting  position  (fig.  2) 

PSICCO  ■  "  0°  ■  eccentricity  angle  of  pallet 

PHIO  ■  139°  -  escape  wheel  starting  angle  of  initial  coupled  motion 

PHICUTD  ■  1485°  «  cumulative  escape  wheel  angle  obtained  from  product  of 
total  engaged  rotor  rotation  and  gear  ratio.  The 
total  rotor  ro:.atlon  for  the  M739  fuze  is  46.41°, 
while  the  gear  ratio  is  32.  Thus,  PHICUTD  »  46.41  x 
32  -  1485° 

MU  -  p  -  0.10  -  coefficient  of  friction  of  gear  train  (pivots  and  tooth- 
to-tooth  contacts)  and  escape  wheel  pivot 

MU1  ■  -  0.10  ■  coefficient  of  friction  of  pallet-escape  wheel  inter¬ 

face  and  pallet  pivot 


Gear  Parameters 


PSUBD1  -  P^i  -  80  -  diametral  pitch  of  mesh  no.  1  (rotor  and  pinion 

no.  2) 

PSUBD2  -  P42  “  100  "  diametral  pitch  of  mesh  no.  2  (gear  no.  2  and 

escape  wheel  pinion) 

NG1  -  Ng1  ■  64  «  number  of  teeth  of  rotor  (full  gear  no.  1) 

NG2  »  Nq2  ■  36  ■  number  of  teeth  of  gear  no.  2 

NP2  “  Np2  ■  9  ■  number  of  teeth  of  pinion  no*  2 

NP3  >  Npj  ■  8  ■  number  of  teeth  of  pinion  no.  3  (escape  wheel  pinion) 
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CAPRF1  ■  Rpj  -  0*41214  in*  (10*468  urn)  ■  pitch  radius  of  gear  no.  1 

(rotor) 

CAPRP2  ■  Rp2  ■  0*19039  in*  (4.836  am)  ■  pitch  radius  of  gaar  no*  2 

RP2  ■  rp2  ■  0.05796  in.  (1*472  am)  ■  pitch  radius  of  pinion  no.  2 

RP3  ■  rp3  ■  0.04231  in.  (1.075  oa)  ■  pitch  radius  of  pinion  no.  3 

(escape  wheel  pinion) 

THETA1  ■  -  24.215°  «  pressure  angle  of  mesh  no.  1 

THETA2  ■  62  ■  27.326°  ■»  pressure  angle  of  mesh  no.  2 

R1  0.250  in.  (6.350  nm)  ■  distance  of  rotor  pivot  from  spin  axis 

R2  “*^2  "  0*317  in.  (8.052  mm)  ■  distance  of  pivot  of  gear  and  pinion 

set  no.  2  from  spin  axis 

R3  "^3  ■  0.309  in.  (7.849  mm)  ■  distance  of  pivot  of  escape  wheel  from 

spin  axis 

r4  -<R^  m  0.304  in.  (7.722  mm)  ■  distance  of  pivot  of  pallet  from  spin 

axis 

RH01  *■  Pj  ■  0.03075  in.  (0.781  mm)  ■  pivot  radius  of  rotor 

RH02  -  p2  ■  0.015  in.  (0.381  mm)  •  pivot  radius  of  gear  and  pinion  no.  2 

RH03  ■  P3  ■  0.015  in.  (0.381  mm)  ■  pivot  radius  of  escape  wheel 

CAP RBI  ■  R^j  ■  0.37588  in.  (9.547  mm)  ■  base  radius  of  gear  no.  1 

CAPRB2  ■  Rjj2  ■  0.16915  in.  (4.296  mm)  ■  base  radius  of  gear  no.  2 

RB2  ■  rjj2  "  0.05286  in.  (1.343  mm)  ■  base  radius  of  pinion  no.  2 

RB3  ■  rjj^  ■  0.03759  in.  (0.955  mm)  ■  base  radius  of  escape  wheel  pinion 

CAPR01  ■  Rqj  ■  0.41425  in.  (10.522  mm)  ■  outside  radius  of  gear  no.  1 

CAPR02  «■  Rq2  ■  0.19404  in.  (4.929  mm)  *  outside  radiuB  of  gear  no.  2 

R02  ■  Tq2  "  0.07670  in.  (1.948  mm)  ■  outside  radius  of  pinion  no.  2 

R03  **  rQj  ■  0.05580  in.  (1.417  mm)  ■  outside  radius  of  escape  wheel 

pinion 

J1  ■  Jj  ■  0  ■  initialization  parameter  for  mesh  no.  1  [The  zero  value 

corresponds  to  earliest  possible  contact  of  mesh  (ref  3).] 

J2  "  J2  ■  0  -  initialization  parameter  for  mesh  no.  2 


Output  Data 


The  date  blocks  following  the  Input  data  represent  the  results  of  various 
computations • 


Fuse  Geometry 


The  angles  BETA1D  -  to  3ETA3D  -  0  and  GAMMA2D  -  y2  Co  GAMMA4D  •  y^ 
are  prlntod  for  checking  purposes* 

Coupled  Motion 

The  first  coupled  motion  output  refers  to  the  entrance  side  of  the 
verge*  For  each  time  T  of  the  coupled  motion,  the  following  variables  are 
computed : 


PUZ  ■  <J>  ■  Instantaneous  escape  wheel  angle  (deg) 

PHIDOT  ■  l  *  escape  wheel  angular  velocity  (rad/sec) 

G  ■  g  ■  pallet  -  escape  wheel  contact  position  (in.)  (equation  C15  of 
appendix  C) 

PSID  *  ■  pallet  angle  (deg) 

PSIDOT  -  ij  ■  pallet  angular  velocity  (rad/sec) 

PHITOT  m  <p,j,  "  cumulative  escape  wheel  angle  (deg) 

F23  *  F23  **  normal  contact  force  of  gear  no.  2  on  pinion  no.  3  (lb) 

F12  -  F12  -  normal  contact  force  of  gear  no.  1  on  pinion  no.  2  (lb) 

PN  »  Pn  «  normal  contact  force  between  escape  wheel  and  pallet  (lb), 
computed  according  to  equation  B185  In  appendix  B 

PNPSI  -  Pn  -  normal  contact  force  between  escape  wheel  and  pallet  (lb), 
computed  according  to  equation  B186  In  appendix  B  (serves 
for  checking) 

*•  2 
DPHI2  »  <f>  ■  escape  wheel  angular  acceleration  (rad/sec  ),  Runge-Kutta 

output 


Free  Motion 


Tha  first  free  notion  on  the  exit  aide  follow*  the  coupled  notion  on  the 
entrance  elde  of  the  verge*  For  each  tine  T  of  the  free  motion,  the  following 
variable*  are  evaluated: 

PHI  ■  $  ■  instantaneous  escape  wheel  angle  (deg) 

PHIDOT  -  $  -  escape  wheel  angular  velocity  (rad/sec) 

PSI  ■  t|»  ■  pallet  angle  (deg) 

PSIDOT  ■  ■  pallet  angular  velocity  (rad/sec) 

PHITOT  «  $T  -  cumulative  escape  wheel  angle  (deg) 

FF12  -  FF12  "  normal  contact  force  of  gear  no*  1  on  pinion  no.  2  for 
free  motion  (lb) 

FF23  ■  Fp23  ■  normal  contact  force  of  gear  no.  2  on  escape  wheel  pinion 
for  free  motion  (lb) 

Impact 


The  first  exit  impact  follows  the  first  exit  free  motion.  Just  preced¬ 
ing  the  "IMPACT"  label,  the  program  prints  the  values  of  VP  -  V~Ni  and  VS  *  VgN! 
which  stand  for  the  pre-impact  velocity  components,  normal  to  the  verge  face,  of 
both  the  pallet  and  escape  wheel  contact  points  (equations  D13  and  D15).  Subse¬ 
quent  to  the  "IMPACT"  label,  the  following  variables  are  evaluated: 


PHI  -  <J>  -  Instantaneous  escape  wheel  angle  (deg),  same  as  before  impact 
PHIDOT  -  |  ■  post-impact  escape  wheel  angular  velocity  (rad/sec) 


PSI  •*  \p  •>  pallet  angle  (deg),  same  as  before  Impact 

PSIDOT  «  ^  ■  post-impact  pallet  angular  velocity  (rad/sec) 

PHITOT  -  ■  cumulative  escape  wheel  angle  (deg),  same  as  before  impact 

VP  *  VTNf  ■  post-impact  normal  velocity  component  of  pallet  at  contact 
point  (equation  D1S) 


VS  -  Vgtyf  •  post-impact  normal  velocity  component  of  escape  wheel  tooth 
at  contact  point  (equation  D13) 


Note  that  in  the  present  program,  the  post-impact  VP  is  equal  to  VS  since  the 
coefficient  of  restitution  is  zero. 
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"Number  of  Turns -to-Arm"  and  Maximum  Contact  Forces 


The  "number  of  turns-to-arm"  at  30,000  rpm  Is  obtained  with  the  help  of 
that  time  T^g^  which  corresponds  to  the  escape  wheel  angle  PHICUTD  ■  1485°. 
Thus , 


number  of  turns-to-arm  -  — - gg-  ■  x  0.05156  ■  25.8 

The  maximum  non-impact  contact  forces  for  the  total  cycle,  for  both 
coupled  and  free  motion,  are  listed  at  the  end  of  the  output. 


CONCLUSIONS 


While  It  was  not  the  purpose  of  this  Investigation  to  undertake  a  parametric 
study  of  the  mechanism  for  which  the  program  was  written,  the  program  was  suffi¬ 
ciently  tested  to  confirm  that  such  a  study  is  possible.  It  may  include  varia¬ 
tions  in  masses  and  moments  of  inertia  of  all  components;  variations  in  the  loca¬ 
tions  of  the  centers  of  mass  of  the  verge  and  the  rotor;  and  variations  of  gear, 
escapement,  and  fuze  geometries;  as  well  as  various  friction  and  coefficient  of 
restitution  conditions. 

Previous  high-speed  motion  picture  observations  of  pin  pallet  escapements 
showed  that  the  impacts  were  essentially  inelastic  and  that,  therefore,  a  zero 
coefficient  of  restitution  was  justified.  Similar  observations  made  on  the 
detached  lever  escapement  of  the  M577  fuze  timer  confirmed  this. 

The  present  work  reports  only  on  a  single  test  run  using  the  M739  fuze  S&A 
data  with  a  system  coefficient  of  friction  of  0.1.  This  is  assumed  to  be  repre¬ 
sentative  of  actual  test  conditions  since  (1)  previous  simulations  of  pin  pallet 
escapements  showed  that  the  range  of  actual  experimental  results  may  be  repro¬ 
duced  with  coefficients  of  friction  between  0.1  and  0.2  and  (2)  a  special  lubri¬ 
cant  is  used  in  conjunction  with  the  M739  fuze.  This  choice  of  coefficient  of 
friction  is  borne  out  by  the  good  agreement  with  experimental  results. 

Noteworthy  is  the  fact  that  the  mechanism  armed  in  approximately  22  turns 
with  a  coefficient  of  friction  of  0.05  and  in  approximately  31  turns  with  a  coef¬ 
ficient  of  friction  of  0.15. 
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APPENDIX  B 

DYNAMICS  OF  ROTOR  DRIVEN  S&A  MECHANISM  WITH  A  TWO-PASS 
INVOLUTE  GEAR  TRAIN  AND  A  VERGE  (PLATE  PALLET)  RUNAWAY  ESCAPEMENT 
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Description  of  Systems  and  Outline  of  Derivations 


This  appendix  gives  derivations  for  the  complete  mathematical  model  of  an 
S&A  mechanlsmt  consisting  of  a  spin  driven  rotor,  a  two  pass,  step-up  involute 
gear  train  and  a  verge-type  runaway  escapement.  Figures  B-l  and  B-2  show  the  two 
types  of  configurations  to  which  this  model  may  be  adapted. 

This  work  uses  geometric  and  kinematic  equations  for  the  verge  escapement 
developed  in  appendixes  C  through  F.  In  addition,  the  dynamics  of  fuse  gear 
trains  is  based  on  the  work  reported  in  reference  B-l.  As  in  previous  efforts  of 
Lowen  and  Tapper  on  the  pin  pallet  escapement  (refs  B-l  and  B-2),  the  following 
three  motion  regimes  of  the  mechanism  are  considered: 


Coupled  Motion 


The  tip  of  the  escape  wheel  is  in  constant  contact  with  either  the 
entrance  or  the  exit  working  surface  of  the  verge  while  it  is  driven  by  the  rotor 
(gear  no.  1)  through  the  gear  and  pinion  set  no.  2.  The  entrance-  as  well  as  the 
exit-coupled  motion  differential  equations,  which  differ  from  each  other  to  some 
extent,  are  expressed  in  the  escape  wheel  variable  <p  and  are  obtained  by  combin¬ 
ing  the  solutions  to  the  Newtonian  force  and  moment  expressions  for  the  individ¬ 
ual  mechanism  components. 


Free  Motion 


The  pallet  and  the  escape  wheel,  gear  train,  rotor  systems  move  Indepen¬ 
dently  of  each  other  in  this  phase  of  the  motion.  Thus,  there  results  one  dif¬ 
ferential  equation  for  the  pallet  in  the  variable  if),  and  another  one  for  the 
combined  system  in  the  escape  wheel  variable  $. 


Impact 


The  formulation  of  the  Impact  regime  is  similar  to  that  shown  In  refer¬ 
ence  B-I  for  the  pin  pallet  escapement.  Again,  the  moment  of  Inertia  of  the 
escape  wheel  and  pinion  no.  3  also  contains  the  referred  mass  properties  of  the 
rotor  and  gear  and  pinion  set  no.  2.  This  impact  simulation  is  based  on  the 
classical  angular  impulse  -  momentum  model,  where  a  coefficient  of  restitution  is 
used  to  account  for  the  energy  losses.  Since  it  is  assumed  that  the  effect  of 
the  impact  force  between  the  escape  wheel  and  the  pallet  is  much  greater  than  the 
effects  due  to  the  driving  torque  of  the  rotor  and  the  various  retarding  torques 
due  to  friction,  the  latter  are  disregarded. 
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Figure  B  1.  Rotor  driven  S&A  device,  configuration  no 


»?:! 


Treatment  of  Friction  Forces 


The  Influence  of  friction  forces  is  considered  both  in  the  coupled  and  the 
free  motion  regimes.  There  is  friction  at  the  Interface  between  the  escape  wheel 
and  the  pallet  during  coupled  motion,  and  there  is  friction  between  the  teeth  of 
the  gears  at  all  of  the  pivots  during  both  of  these  types  of  motions.  As  in 
reference  B-l,  the  individual  pivot  friction  torques  are  obtained  by  the  alge¬ 
braic  addition  of  the  two  friction  moments  due  to  the  x  and  y  components  of  the 
normal  bearing  forces,  rather  than  by  the  direct  use  of  the  resulting  normal 
forces.  This  conservative  approach  to  friction  is  necessary  to  avoid  the  diffi¬ 
culties  which  the  presence  of  a  square  root  term  Introduces  into  the  solutions  of 
the  various  differential  equations. 


Appendix  E  of  reference  B-l  proves  that  a  change  in  sign  of  the  coefficient 
of  friction,  both  in  the  coupled  and  the  free  motion  simulations  of  the  gear 
train,  including  the  rotor  and  the  escape  wheel  (with  the  exception  of  the 
escape  wheel-pallet  Interface),  is  sufficient  to  account  for  the  changed  direc¬ 
tions  of  the  friction  forces  encountered  during  a  possible  motion  reversal  of  the 
gear  train  following  an  impact. 


Geometry  of  Fuze  Body  Configurations 


The  following  section  contains  derivations  of  expressions  for  the  various 
fuze  body  angles  associated  with  the  pivot  holes  of  both  configurations  of  the 
two  pass,  step-up  gear  mechanisms.  In  addition,  relationships  between  the  unit 
vectors  of  the  body-fixed  X-Y  and  x'-y*  systems  (figs.  B-l  and  B-2)  are  given  to 
describe  certain  escapement  forces  in  terms  of  the  primed  coordinate  system. 
(Note  that  the  x'-axls  is  oriented  along  the  escapement  center-line  0  ~0S.) 
Finally  a  slgnum  function  is  introduced  so  that  common  programming  expressions 
for  both  configurations  can  be  devised. 


Fuze  Body  Configuration  No.  1 


The  following  relationships  for  the  angles  Y^,  <5^,  and  8^  for  configura¬ 
tion  no.  1  are  indicated  by  figure  B-3: 
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Further , 


-1  &  +(il  ~  (Rp2  +  r»3>‘ 


Y3  -  CO*’1  [ 


Finally, 


Y3  ■  y2  +  Y3 


The  angle  y,  1*  given  by: 


Then, 


•  -1  ,*!  +fi4"  *2, 

Y«  ■  “>•  [ — 3T^ - 1 


y4*y3+y4 


Angles  6^.  The  angle*  6^  and  63,  a*  shown  In  figure  B-3,  become: 
-1  <Rnl  +  ro2)2  +  ’fl2, 

s* ' co<  ft'"#'  1 


-1  <Rn9  +  r„^>2  +fl9 

V~‘ 


The  angle  6^  must  take  the  pivot-to-pivot  distance  a  of  the  escapement  Into 
account •  The ref ore , 

.  a2 

64  -  cos  [ - 53J- - -]  (B8) 

Angles  8^.  The  angles  8^  and  6?  become  respectively: 

0A  -  tt  -  «2  (B9) 


S2  "  Y2  +  11  "  *3 


(BIO) 
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(BU) 


The  angle  ft.  is  found  with  the  help  of  angle  6  of  equation  B8»  Then, 

3  ^ 

1*3  "  Y3  +  *  -  \ 

Figure  B-3  shows  the  angle  y’  between  the  positive  x'-axls  and  the  unit  vector 
i»4.  It  Is  given  by:  p 

Y*  -  ir  -  *  (B12) 

P 

where  e  Is  obtained  with  the  help  of  equations  B4  and  B8: 

e  -  it  -  64  -  yJ,  (B13) 

and  therefore: 

Y'  -  64  +  Yj  (B14) 

The  unit  vector  n,  is  expressed  in  terns  of  the  primed  coordinate  system  as 
follows: 

n.  ■  cos  Y +  sin  y'J'  (B15) 

4  P  P 

Further,  the  unit  vectors  7'  and  3',  when  expressed  in  the  X-Y  system,  become: 

1’  -  cos  (ft3  -  ir)l  +  sin  (ft3~  ir)J 
or 

1'  -  -cosftjl  -  sinft3j  (B16) 

end 

J*  -  k*  x  I*  -  sinftjl  -  co8ft3J  (B17) 

Fuze  Body  Configuration  No.  2 

The  angles  y.,  5  ,  and  ft.  of  configuration  no.  2  are  defined  in  figure 
B-4.  11  1 

Angles  y^*  Since  these  angles  are  defined  in  the  clockwise  direction 
with  respect  to  the  body-fixed  X-axis,  their  values  must  be  negative. 


jr  —  w  *  •  *%.',  •v*Tj 


Then, 


Y2  -  -CO.  I  - *  J 


(B18) 


Anglo  y!  1.  givon  by  oquotlon  B2  and  with  that  tho  following  oxpreaslon  1. 
obtained1* for  y^t 


Y3  -  Y2  -  Yj 


(B19) 


Tho 


nogativo  sign  for  y3  la  necoaaary  In  ordor  to  mako  Y3  negative • 


Tho  angle  yI  1.  given  by  equation  B4.  It  allows  the  following 
expression  for  angle  Y^t 


Y4  -  Y3  -  Yj 


(B20) 


Again,  y^  must  have  a  negative  sign  to  make  Y^  negative. 


Angles  <5^.  The  angles  6_ »  <s3»  «®d  6^  can  be  taken  directly  from  equa¬ 
tions  B6,  BY,  and  B8,  respectively. 


Angles  0^.  Angle  0^  becomes : 


Further 


Finally, 


*1  "  "  +  fi2 


02  "  Y2  +  *  +  «3 


(B21) 


(B22) 


% 


S3  "  Yj  +  ■  +  «4 


(B23) 


The  angle  Y’  between  the  positive  x’-axls  and  the  unit  vector  ii ’  is 
shown  in  figure  B-4.  This  angle  is  given  by:  4 


where 


Therefore 


Y*  ■  ir  +  e 

P 


e  "  *  -  64  "  Yi 


(B24) 


(B25) 


V' 

i 


f 

r 


Yp  •  2it  -  6^  -  y^ 


(B26) 
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-AL- 


Tht  unit  Victor  n.  can  again  ba  expressed  in  tarns  of  tha  primed  coordinate  ays- 

tan  as  * 


nA  «  cos  Yii'  +  sin  y'J' 

“  P  D 


(B27) 


follows i 


Tha  unit  vectors  I’  and  J '  are  expressed  in  tha  X-Y  system 


i*  *  -  cos  0.1  -  sin  0.J 


(B28) 


J'  -  sin  63!  -  cos  83J 

Common  Computational  Expressions  for  Both  Configurations 


(B29) 


_  .  f*nd  common  programming  axprasslons  for  all  angles  of  both  configura¬ 
tions  (with  tha  exception  of  tha  angles  y*  of  equations  B14  and  B26)  the  follow¬ 
ing  slgnum  function  Sg  la  Introduced!  p 


eg  ■  +1  for  configuration  no*  1 
Sg  ■  -1  for  configuration  no*  2 
This  leads  to  the  following  expressions! 

According  to  equations  B1  and  B18! 


(B30) 

(B31) 


-1  -  OL ,  +  r  y 


According  to  equations  B3  and  B19: 


Y3  -  y2  +  s6  y3 


(B34) 


According  to  equation  B4: 


<R?  +fl2  -  M2 


According  to  equations  B5  and  B20» 


Y4  -Y3  +  ,6  Y4 


(B36) 


According  to  equations  B6,  B7,  and  B8,  raapactlvaly 


52  *  cot 


6,  »  cos”A  [ 


-i  r<yy2+fl?-< 

I  aV'*pi 

-1  p2  _  rp3j  *  2  3i 


+  rp3) 


-1 

co*  t — 23 r — 


] 


(B37) 


(B38) 


(B39) 


Dynamics  of  Pallet  and  Escape  Wheel  in  Coupled  Motion 


This  section  shows  derivations  for  both  entrance-  and  exit-coupled  motion 
differential  equations.  While  the  resulting  combined  differential  equations  are 
similar,  they  differ  in  a  number  of  terms  because  of  the  inherent  differences  of 
the  nt  -  n^  coordinate  systems  for  the  entrance  and  exit  sides.  While  these 
derivations  are  given  for  configuration  no.  1,  they  are  also  applicable  to  con¬ 
figuration  no.  2  with  the  help  of  the  slgnum  function  8g  introduced  in  the  pre¬ 
vious  section. 

The  analysis  of  the  pallet  is  performed  in  the  x'  -  y'  system  for  conve¬ 
nience,  while  that  of  the  escape  wheel  uses  the  X  -  Y  system.  An  appropriate 
coordinate  transformation  makes  it  possible  to  combine  the  systems. 

The  coefficient  of  friction  u.  is  used  at  the  escape  wheel-pallet  interface 
as  well  as  at  the  pallet  pivot,  while  the  coefficient  of  friction  u  serves  to 
express  the  friction  forces  at  the  escape  wheel  pivot  as  well  as  at  the  escape 
wheel-pinion  no.  2  tooth  interface,  (v  is  also  used  in  all  other  gear  meshes  and 
gear  pivots.) 


Slgnum  Functions  and  Acceleration  of  Center  of  Mass  of  Pallet 


To  define  the  direction  of  the  friction  force  of  the  escape  wheel  point 
S  on  the  coincident  pallet  point  T,  the  slgnum  function  s4  is  introduced.  Thus, 


s4 


VS/T 

Ivs/tI 


(B40) 
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whir*  Vg  v«  la  glvan  by  aquation  C2S  of  appendix  C.  Tha  aignua  function  a.  aarvaa 
for  the  'determination  of  tha  dlractlon  of  tha  friction  forcaa  acting  on  tha  pal- 
lat  pivot  *  Than 


a,  -  -4-  (B41) 

5  Hi 

With  tha  conatant  apin  w  of  tha  fuaa  body,  tha  accalaratlon  X  of  tha  cantar  of 
aaaa  of  tha  pallat  la  glvan  by  (figs.  B-l  and  B-5):  cp 

Acp  ■  -w2^n4  -  <w  +  i)2  rCpIcoa(i|»  +  q»c>i*  +  ain<*  +  i|fc)I'J 

+  J  rcp  [-ain(i|»  +  4>c)i'  +  coa(<|>  +  t|»c > J ’ ]  ^B*2^ 


Dynamics  of  Pallat  and  Bscape  Whaal  for  Bntranca-Coupled  Motion 


A  datailad  fraa  body  diagram  of  tha  pallat  for  entrance-coupled  motion 
la  shown  In  flgurs  B-5.  The  free  body  diagram  of  tha  escape  wheel  is  simplified 
and  presented  for  orientation  only*  A  complete  free  body  diagram  of  the  escape 
wheel  which  includes  all  pivot  and  contact  forces  is  given  by  figure  B-6. 


Force  Equation  of  Pallet .  Substitution  of  all  forces  into  Newton's  Law, 
together  with  the  acceleration  of  the  center  of  mass  according  to  equation  B42 
furnishes: 


V,  +  V4pn"t  -  'ip  l'  -  VsV'  +  V  +  V/y 


yp“ 


xp- 


-  {-u)^«4n4  -  (w  +  $)2rcp  lcos(i|/  +  <|>c)i'  +  sin(tp  +  ii»c)j'] 

+  Y  rcp  I-sin(ij/  +  ij»c)i»  +  cos(t|>  +  ^C)J'J> 


(B43) 


(For  nomenclature,  consult  figures  B-l  and  B-5,  as  well  as  reference  B-l.) 


Moment  Equation  of  Pallet*  The  moment  equation  of  the  pallet  must  be 
written  wTtK  the  understanding  that  pivot  point  0  is  an  accelerated  point. 
Thus ,  p 


where 


rop  x  mp  rcp 


+  H 


op 


(B44) 


■  sum  of  all  external  moments  with  respect  to  pivot  Op 
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absolute  acceleration  of  point  0p 


H  *  tlaa  rate  of  change  of  angular  momentum  of  the  pallet  with 
respect  to  point  0p. 

Since  the  spin  angular  velocity  u>  of  the  fuse  body  is  constant!  equation  B44 
takas  the  following  form: 


M  -  -(-u^R.n-)  x  n»  r  Icoa<ij»  +  *  )i*  +  sin(i|>  +  ib  )J*]  +  (B45) 

op  p  cp  c  c  p 


If  n^  is  substituted  according  to  equation  B27,  the  above  becomes  in  scalar 
terms: 


M  -If-  a i  r  fl.w2sin(Yfi 

Op  p  p  Cp  4  p  C 


(B46) 


To  find  the  value  of  f^p  consider  figure  B~5. 


"op  ■  Di"t  *  Vn  +  Cl”n  x  WA  -  »pul*5(,xp  +  V*  W7) 


or,  in  scalar  terms 


M  -  DIP 
op  In 


*  Ci“l*«pn  '  pl“l*5<Fxp  +  Pyp) 


(B48) 


When  equation  B48  is  substituted  into  equation  B46,  the  final  moment  expression 
for  the  pallet  is  obtained: 


P  ID’  -  C'p.s,]  -  p  p.s-CF  +  F  ) 
n  1  1  11  4J  p^l  5  xp  yp 


I  ip  -  m  fi  r  u)Zsin(Y,!  -♦-*_) 
p  p  4  cp  p  c 


(B49) 


As  originally  discussed  in  an  ARRADCOM  report  (ref  B-3),  Fxp  and  F^.p  represent 
conservatively  evaluated  pivot  force  components  which  assure  that  the  pivot  fric¬ 
tion  moments  are  opposed  to  the  rotation  at  all  times. 


The  pivot  friction  components  F^  and  F^p  are  first  obtained  from 
the  following  component  expressions  of  equation  B43.  Subsequently,  they  are 

mi  mi 

transformed  to  become  F  and  F  ,  respectively. 

xp  yp 


-Pnsin(*  +  a)  +  u1s4Pncos(^  +  a)  +  F^p  -  WjS^ 


2  t  *  2  " 

-  mp(-u>  fl4co8Yp  -  (<j  +  ip)  rcpcoe(f  +  d»c)  -  <j>rcpsin(i|/  +  <J»c)] 


(B50) 
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PnC0t(iJ>  +  ct)  +  +  a)  +  wi,5Fip  +  Fyp 


(B51 ) 


Pp  (-w2if^sinY^  -(<*>  +  ^)2rcp*in(t  +  \|»c)  +  t|»rcpcos(\|>  +  ^C)J 


Simultaneous  solution  of  ths  sbovt  furnishes t 


*„  -  Vi +  “%  *  ai5r  S  *  JS  *  *\ 


(B52) 


S  "  prA  +  A*  *  2  K  *  A-  *  iAft 

Xp  n  5  6  J3J*  7  7  8 


(B53) 


where 


+  s j)sln(ij/  +  a)  -  <1  -  s^s5upcos(^  +  a) 

Ai "  I - TT~2 - 1 

1  +  Mj 


(B54) 


m  («,(slnY'  -  u  s  cosy')  +  r  (sln(t|»  +  t|»  )  -  u.s.cosC^  +  t|»  ))] 
P  *t _  P  15 _ p  cp _ C  1  5  c 


1  +  y 


(B55) 


o  r  [sin(<Ji  +  *  )  -  y  s  cos(i|>  +  *  )) 
P  cp _ c  15  c 


1  +  y 


(B56) 


m«r*.Jcos^  +  O  +  ♦  )J 

|  P  CP  c  1  5  C  | 

1  +  y? 


(B57) 


i  yl^84  +  +  u)  +  <1  -  s485ypsin(\f»  +  a) 

I  2  '  I 


CB58) 


l  +  y; 


mn^“flA^C0*Yn  +  ui8s8lny«)  “  r  (cosC  1|»  +  ♦  )  +  y.s  Jln(t  +  ♦  ))J 

a6  -  |  J— i - 2 - U - E - SfZ - £ - U - £_  |  <M9) 

1  +  yf 


(B60) 


mr[cos(ij)  +  *  )  +  +  *  )] 

A,  .  I  -i-SB - c  ■  •»  1  5 - -  | 


*8  ' 


1  +  u 


1 


a r  lein(i|i  +  *  )  -  u.a-coaC^  +  i|>  )J 

P  CP  C  X  J  c 

1  +  M? 


(B61) 
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The  factor  u>  /|u>|  la  Introduced  in  place  of  u  to  assure  that  the 
quantity  Is  positive  regardless  of  direction  of  spin.  (This  too,  assures  that; 
the  friction  moments  oppose  rotation*  The  driving  moment  is  proportional  to  u 
and,  therefore,  independent  of  spin  direction*)  To  make  the  final  decision  con¬ 
cerning  the  signe  in  equations  B52  and  B53,  these  forces  are  now  substituted  into 
moment  equation  B49,  and  the  influence  of  the  direction  of  rotation  on  each  of 
the  resulting  moment  computations  is  explored: 

Pn*Di  “  C1M184  "  PpUl85^Al  +  A5^  *  ppWl85u2^A2  +  A6^ 

*  2pppl85R^<A3  +  V  *  Vl85*2(A3  +  V  ±  ppMlV(A4  +  V  (B62) 

-  y  -  mprcp<R5a,28ln(Yp  ■  *  -  v 


With  s5  positive  for  positive  rotation,  and  vice  versa,  and  with 
all  other  parameters  positive  at  all  times,  the  following  moment  components  of 
equation  B62  must  have  negative  signs  during  positive  rotation: 


-VpVs^i  +  V 

2 

-PpMjSjO)  (A2  +  Afi) 


(B63) 

(B64) 


V 


(A65) 


the  sign  of 
becomes : 


The  sign  of  the  term  containing  ^  is  negative  and  is  controlled  by 
Therefore,  the  slgnum  operator  a<-  is  omitted  and  the  term 


(B66) 


The  choice  of  sign  for  the  friction  moment  term  in  equation  B62, 

*• 

which  is  proportional  to  the  pallet  angular  acceleration  ij>,  is  discussed  in 
detail  in  appendix  F  of  reference  B-l*  This  work  leads  to  the  computational 
rules  of  equations  B73  and  B74  below  which  deal  with  the  sign  in  the  expressions 
for  the  effective  moment  of  inertia  IpR*  (The  slgnum  function  s^  has  now  been 
omitted.) 
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With  these  considerations,  equation  B62B~*  now  becomes 


Vis  -  "2a19  '  *R**20  -  i  *21  (B68) 

■  *»♦  '  ”prcpfl4“2aln<Yi  '  *  -  V 

where®"^ 

*18  ■  Bi  '  Cl“l*4  “  Vl*5(*l  +  *5)  <B69) 

*19  *  pplll*5^*2  +  *6^  (870) 

*20  '  Vl<A3  +  V  <B71) 

*21  '  pp“l*5<A3  +  *7>  <B72) 

•  ## 

IpR  ■  I  +  A22*  when  ij>  and  $  have  Identical  signs  (B73) 

•  ** 

IpR  -  Ip  -  A22*  when  t|>  and  <j>  have  opposite  signs0  (B74) 

A22  "  ppyl(A4  +  V  (B75) 


Equation  B68  is  now  used  to  find  an  expression  for  the  contact 
force  P  .  This  expression  will  later  make  it  possible  to  establish  a  single 
differential  equation  for  the  escapement  in  coupled  motion.  Thus, 


B~1 

Note  that  there  is  no  equation  B67. 

The  parameters  An  are  not  sequential  at  this  point.  Ag  to  A.y  may  be  found 
in  equations  B90  to  B95  and  B97  to  B99  in  conjunction  with  the  work  on  the 
escape  wheel* 

Care  must  be  taken  that  I  -  A22  does  not  become  negative*  If  this  occurs, 
IpR  must  be  set  equal  to  zero* 
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It  is  now  necessary  to  write  equation  B76  in  terms  of  the  angular  velocity  $  and 

•• 

the  angular  acceleration  <j>  of  the  escape  wheel  so  that  it  may  later  be  equated  to 
a  similar  expression  for  the  escape  wheel* 

Thus,  with  the  help  of  equations  Cl 9  and  C26  of  appendix  C,  the 
contact  force  becomes: 

Pn  '  A^  1W»  +  (A21u2  +  W>*2  +  *R  A20"»  (B77) 

+  d2A1}  -  »prc(tf4w2«ln<rp 


Escape  Wheal  Analysis*  A  detailed  free  body  diagram  of  the  escape  wheel 
and  pinion  no.  3  in  entrance-coupled  motion  is  shown  in  figure  B-6.  The  pivot 
forces  P  j  and  F  j  as  well  as  the  contact  force  F23  and  the  force  T3  are  now 
defined  in  the  general  X-Y  system.  This  makes  it  necessary  to  transform  the  unit 
vectors  nc  and  ^  from  the  x'  -  y*  system  to  the  X-Y  system.  Since 


or 

i '  -  cos ( $3  -  ir)  1  +  sin(P3  -  ir)J 

(B78) 

and 

I*  a  -cosfjjl  -  sinBjj, 

(B79) 

J'  -  sinful  -  cosB3J 

(B80) 

the  unit  vectors  cL  and  can  now  be  determined  with  the  help  of  equations  C5 
and  C6  of  appendix  C.  Thus, 


-cos(i|>  +  a  +  $3)i  -  sin(\|;  +  a  +  &3)J 

(B81) 

sin(ip  +  a  +  S3)i  -  cos(\J>  +  a  +  S3)J 

(B82) 

Parallel  work  is  given  in  appendix  B  of  reference  B-l. 


Force  Equation  of  Escape  Wheel.  Force  equilibrium  is  obtained  with  the 
help  of  the  following  expression,  based  on  figure  B-6: 


■pnnn”ul*4Pnnt  +  F23n23  +  M*2F23nN23  +  T3 


+  Fx3i  +  ,:Fw^i  +  uF„,j  -  F„,j 


y3 


x3J 


y3J 


(B83a) 
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If 


h 


ik;  .'t 

[i 


The  aignun  function  s«  la  defined  In  in  reference  B-3.  Further,  according  to 
reference  B-l,  the  D'Alembert  force  T3  la  given  by 

T3  "  m3fl3M  (C0*Y3i  +  «inY3J)  (B83b) 


Moment  Equation  of  Bacape  Wheal*  Again,  baaed  on  the  nomenclature  of 
figure  B-6,  the  moment  equation  of  the  eacape  wheel  la  given  by: 


AlV("Pn,%  +  Vn  * 

-U|)3<Px3  +  Py3)E  +  rb3F23E  '  “*2<d2  '  *2>F23E  '  VE 


(B84) 


This  becomes  In  scalar  form: 

‘Pn<\  -  Vl‘(>  "  »3<F*J  +  V  +  rb3P23 
-u,2(d2  '  *2)P23  '  V  (,8!> 

*w  m 

The  pivot  force  components  Fx3  and  Fy3,  which  eventually  become  and  Fy3>  are 

now  determined  with  the  help  of  the  following  scalar  component  expressions  of 
equation  B83: 


Pnsin(i|i  +  a  +  03)  +  ^s^cos^  +  a  +  8-j)  +  T3cosy3 
+  F23sin(82  +  02)  +  us2F23cos(82  +  02>  +  Fx3  +  MFy3  -  0  (B86) 


Pnc08(i|>  +  a  +  B3)  +  u1s^Pnsin(<|»  +  a  +  83)  +  T38inY3 
“F23cos(82  +  62)  +  us2F23sin(82  +  02)  +  pFx3  -  Fy3  -  0  (B87) 

Simultaneous  solution  of  the  above  yields  the  following  "conservative"  pivot 
force  components : 


Fy3  "  ^*9  +  F23A10  +  T3AU 


Fx3  "  PnA12  +  F23A13  +  T3AU 


(B88) 


(B89) 
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Combined  Entranced- Coupled  Motion  Differential  Equation .  Equations  B77 
and  B96  ara  now  sat  equal  to  each  other.  Vhia  furnishes  the  following  comblnad 
motion  dlffarantlal  aquation  for  tha  aacapamant  t 

!IiAie  +  +  ‘'pa*!;7  +  +  (,100) 
*  A13Aie'23  -  A16A18T3  "  A17“8[A19  '  ‘  *  *  V1 


Dynamics  of  Pallet  and  Escapa  Wheel  for  Exit -Coupled  Motion 


A,  free  body  diagram  of  tha  pallet  for  exit-coupled  motion,  together  with 
a  simplified  representation  of  tha  escape  wheel  la  shown  in  figure  B-7.  When 
this  figure  Is  compared  with  figure  B-5,  which  depicts  entrance-coupled  motion, 
it  is  noted  that  now 

P„  -  -P_n  (B101) 

n  n  n 

This  sign  change  Is  reflected  both  in  the  force  and  In  the  moment  equation. 


Force  Equations  of  Pallet.  When  the  scalar  force  equations  B50  and  B51 
are  altered  to  account  for  the  above  change,  they  have  the  following  form: 


Pnsin(t  +  a)  +  UjS4Pncos(i|>  +  a)  +  F^p  -  jjjSjF^  (B102) 

2  *2  •• 

-  Sp[-u  <R^cosYp  -  <w  +  \|»)  rcpcos(\|*  +  i|>c)  -  i|/rcpsln(i|j  +  \|>c)] 


— PjjCos  (ty  +  a)  +  UjS^sin^  +  a)  +  ^SjF^  +  F^p  (B103) 

-  ®p I-u)^(R^sinYp  -  (u>  +  |)2rcpsin(t|/  +  +  i|>rcpcos(<|i  +  ipc ) J 


Simultaneous  solution  of  the  above  shows  that  only  the  factors  A.  and  A^  as  given 
by  equations  B54  and  B58  need  be  changed.  These  will  now  be  called  AA,  and  AAy 
Then, 

2 

U,(s,  -  s,  )sin(ij)  +  a)  +  (1  +  p.s.s  ,)cos(i|/  +  a) 

Mi  -  |  -Ui - £ - - - LL3. -  |  (8104) 


and 


AA5 


wl(s5 


-  s ^)cos (i|>  +  a)  -  (1  +  u^s^s j)sin(iJ;  +  a) 


1  +  u 


2 

1 


(B105) 
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Mo—  nt  Equation  of  Pallet*  Tho  sign  change  of  tha  contact  fore*  P,, 
reflects  itself  in  tha  acalar  — nt  aquation  B49  aa  follows: 


“n1"1!  *  Ci“l,4<  *  Vl*3<Fxp  +  V 

p  p  a  cp  p  c 

(Nota  change  of  sign  of  D[») 


(B106) 


Following  tha  aa—  path  aa  for  tha  entrance  contact.  It  can  be 
shown  that,  with  tha  exception  of  tha  factor  Ajg,  aquation  B76  for  the  contact 
force  PQ  re— Ins  unchanged.  Then, 


.  V  *  A21»2  +  2  tSt*20*  +  A,  '  V<a.flX,tl>(lfi  ~  ♦  ~  V 


Iff 


18 


where 


AAia  ■  -IDJ  +  ci»J1»4  +  ppWi*5(AAl  +  M5^ 


(B107) 


(B108) 


Finally,  the  expression  equivalent  to  equation  B77  becomes t 


Pn  ■  *4  f1™0*  +  <A2»"2  +  IP«V,{2  +  2  M  A20Ui 
2  2 

+  0)  A19  -  oprcpfl4u  ain(Yp  -  *  -  ] 


(B109) 


Force  Equations  of  Escape  Wheel*  A  free  body  diagram  of  the  escape 
wheel  in  exit-coupleid  motion  is  shown  in  figure  B-8.  When  the  scalar  force  equa¬ 
tions  B86  and  B87  for  the  escape  wheel  in  entrance-coupled  motion  are  adjusted  to 
account  for  this  condition,  the  following  expressions  are  obtained: 


Pnsin(i|>  +  a  +  +  MjB^P^cosCiJ/  +  a  +  8^)  +  T^cosYj 

+  F23sin(82  +  02)  +  ys2F23cos(02  +  02)  +  Fx3  +  pFy3  -  0  (B110) 


-P  cos(\j>  +  a  +  0.)  +  u.s.P  sin(i 1/  +  a  +  0_)  +  T_8inY- 
n  j  i  4  n  i  5  j 

-F23cos(02  +  02)  +  MS2F23sin(02  +  62>  +  uF^  -  Fy3  -  0  (BUI) 
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Simultaneous  solution  of  ths  above  for  F  «  and  P  j  shows  that  certain  coeffi¬ 
cients  of  these  pivot  forces  will  differ  from  those  given  in  the  solutions  for 
entrance-coupled  motion.  Thus, 


+  F23A10 


+  Vn 


(B112) 


■  P-AA 
n 


12 


+  P23*L3 


+  Vl* 


(B113) 


where  A^q,  A^,  A^,  and  A^  are  given  by  equations  B91,  B92,  B94,  and  B95, 
respectively.  The  new  coefficients  AAg  and  AA^  are  defined  ae  follows: 

(p.s.  -  p)ain(i|>  +  a  +  B,)  -  (1  +  s.pp  )cos(i{/  +  a  +  8  J 

AAq  -  I  -1-5 - 3-5 - 5—1 - 3-  |  (DUO 

*  1  +  u 


-(1  +  pp.s , )sin( ip  +  a  +  B,)  +  (p  -  p.s.JcosC*  +  a  +  8,) 

A*,,  -  | - 1-5 - 3_ - L* - 3_|  (BUS) 

1  +  p 


Moment  Equation  of  Escape  Wheel.  Modification  of  the  escape  wheel 
moment  equation  B&4  for  exit  contact  gives: 

A{  "t  X  Vn  +  51  x  (“Mla4Pn>;it  “  MP3(Fx3  +  Py3)E  (B116) 

•  • 

+  rb3F23^  “  ys2^d2  “  a2^F23^  "  I3^ 

Substitution  of  equations  B112  and  B113  into  the  scalar  form  of  equation  B116 
above  furnishes  the  following  expression  for  Pn: 


Pn  - 


~V  +  F23A15  "  T3A16 


AA 


17 


(B117) 


where  Ajg  and  Ajg  are  given  by  equations  B97  and  B98,  respectively.  The  factor 
AA^j  is  new;  i.e., 


AA17  -  -A[  -  w1s4B'  +  pp3(AA9  +  AA12) 


(B118) 


Combined  Exit-Coupled  Motion  Differential  Equation.  Equations  B109  and 
B 1 1 7  are  now  equated  to  each  other.  This  furniskes  the  combined  exit-coupled 
motion  differential  equation: 
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■  *15*^8*33  '  A16AA18T3  “  A4l7“2<A19  "  Vcp^,ln(YP  '  *  *  *c)J 

(B119) 

Dynamics  of  Rotor  (Goar  No.  1) 

(Applicable  to  both  configurations  with  appropriate  choice  of  Sg.  See  also 
reference  B~l.) 


A  free  body  diagram  of  the  rotor  of  configuration  no.  1,  which  moves  In  a 
counterclockwise  direction,  is  shown  In  figure  B-9.  The  acceleration  of  its  cen¬ 
ter  of  mass  is  given  by  equation  A122  of  appendix  A  of  reference  B~l. 

Since  the  motion  of  the  rotor  must  be  expressed  in  terms  of  the  escape  wheel 
variable  </>,  use  must  be  made  of  the  following  relationships: 


*1  *  N31* 

(B120) 

and 

•* 

*1  ■  *31* 

(B 1 2 1 ) 

where 

„  “P3  «  NP2 

31  Ng2  x  NC1 

(B122) 

The  rotor  angle  <j>lc  +  <p^  Is  expressed  in  terms  of 

V 

the  total  escape  wheel  rotation 

^1C  +  *1  “  *1C  +  N31  ^T 

(B123) 

(For  details, 

see  Program  Description.) 

Newton's 
reference  B-l 

force  equation  may  now  be  written 
.  The  tooth  forces 

in  the  form  of  equation  B63  of 

F21  “  "F12  n12 


(B124) 


and 


Ff 21  "  M81  F12  nN12 


(B125) 


have  Che  appropriate  directions  to  account  for  friction.  (The  signum  function 
is  defined  in  reference  B-l.) 

Thus: 


where 


and 


”F21n12  +  U8lF12nN12 
2  - 

■  nij  {-w  <$^i  -  (w  + 


+  N31*rclI"8in(*lc 


+  FxlA  “  ^yl1  +  FylJ  + 

N31*)2rcl[c08(^c  +  N3lV*  +  8in  (*lc 
+  N2^T)i  +  cos($lc  +  N31d>T)J]> 


+  VVJl 

(B126) 


n12  -  -sin(01  -  0^1  +  cos  ( -  O^J 

(B127) 

“N12  "  “C08^i  “  •“  BinCSj  “  )J 

(B128) 

The  moment  equation  must  be  written  with  respect  to  the  accelerated  point 

M  •» 

Oj.  The  pivot  reactions  Fxi  and  F  j  are  written  such  that  the  associated  fric¬ 
tion  moments  retard  the  counterclockwise  rotation  of  the  rotor*  This  leads  to: 


-F12Rbl  +  M8lalF12  "  upl(Fxl  +  Vyl) 

■  V^'cl*1**!*  +  N3lV  + 


(B129) 


The  forces  F  .  and  F  .  are  obtained  after  the  simultaneous  solution  of  the 
following  component*  expressions  of  equation  B63  of  reference  B-l  for  Fxj  and  F^: 


and 


FX2,la(Sl  "  V  “  “■xF12cob<B1  -  V  +  F*1  -  uFyl 

-  -  (>  +  H3i;)2rclco.Uu  +  N31*t)  -  N3l5rcl.ln(*lc  +  N^)] 

(B130) 


-F12co.(Sj  -  Bj)  -  B.jFjj.lnXSj  -  8j>  +  FyJ  +  uF,, 

-  *,[-(<»  +  »31i)*rcl.l»(4lc  +  K31*t)  +  N3iJrclco.(*lc  +  »3,.T)I 


(B131) 
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*  fn  ih 


i:*"VTVi  *  v'K7r',v'*;B  rr'*' 


m 


Equations  B132  and  B133  are  now  substituted  into  the  moment  equation  B129. 
“F12Rbl  +  M8ialF12  “  PP1[F12(A23  +  A27J  *  “  (A24  +  A28J 


*  2FT  N31^A25  +  A29^  *  ^N3i^2^A25  +  A29^  *  N31**A26  +  A30^ 

’  vVcl8ln(<l>lc  +  N31 V  +  X1N31* 


(B142) 


In  order  to  have  the  pivot  friction  moment  negative,  the  following  terms  on 
the  left  hand  side  of  equation  B142  must  be  negative: 

WP1  F12  (A23  +  A27J 


MPl  a)  <A24  +  A2g) 

“l  (1V)2  <A25  +  V 


2lulpl  KB31J  <A2S  +  A29> 


The  term  up.N.  4>(A  +  A  )  is  treated  in  the  same  manner  as  was  shown  in  connec¬ 
tion  with  equation  *147  or  reference  B-l*  (The  discussion  in  appendix  A  of  ref¬ 
erence  B-l  is  also  of  interest  for  the  determination  of  the  effective  moment  of 
inertia  I^R.) 

The  above  considerations  give  the  moment  equation  the  following  form: 


"F12tRbl  “  U8lal  +  MP1<A23  +  A27)J  "  “  ypl(A24  +  A28) 


2WPJU1  N31<J) 


(A25  +  A2j)  -  iipj(N3j£)  (A23  +  A2j) 


"  v  Vcl8in(*lc  +  N3lV  +  lIi  *  MP1<A26  +  A30),N31* 


(B143) 


With  N»,  <j>  ■  positive,  must  be  absolute.  See  also  discussion  following  equa¬ 
tion  AI47  of  reference  B-l, 


I 


Finally*  the  above  expression  Is  solved  for  F.2,  using  the  effective  moment 
of  Inertia  I^Ri 

-A32  -  A33N3i$  ”  A34^N3i^  “  A35sln<tic  +  “  IlRN31<f> 


(B144) 

where 

A3X  "  Rbl  “  y8lal  +  wpl^A23  +  A275  (B145) 

A32  "  ppiu2<A24  +  A28>  (B146) 

A33  "  2lpK  T^T  <A25  +  A29)  (B147) 

A34  *  MP1  (A25  +  A29>  <B148) 

A35  *  “lVcl“2  (B149) 


I1R  *  I1  +  lp!pl^A26  +  A30^*  wlien  ♦  an<*  ♦  have  the  8ame  signs*  (B150) 

#  •• 

I^R  ■  Ij  -  |u|pj  (A2g  +  A3q)»  when  4  and  <p  have  opposite  signs.  (B151) 

Dynamics  of  Gear  and  Pinion  Set  No.  2  (applicable  to  both  configurations  with 
appropriate  choice  of  8g ) 

The  free  body  diagram  of  gear  and  pinion  set  no.  2  Is  shown  In  figure  B- 
10.  Its  force  equation  Includes  the  D'Alembert  force 

T2  -  mfy2  (B152) 

Thus, 

"F23n23  “  82uF23”n23  +  F12n12  "  y8lF12"NI2  +  Fx2*  +  uFy2*  (B153) 
+  Fy2J  “  wFx2J  +  T2(cosY2i  +  slnY2J)  -  0 
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where 


n23  -  sin(U2  +  e2)i  ■-  cos  (02  +  02)J  (B154) 

and 

"N23  “  coa<&2  +  V*  +  8ln<62  +  02)J  (B155) 

The  unit  vectors  n12  and  nK12  were  given  by  equations  B127  and  3128.  Hence 

*32  "  “F23  "23  (B156) 

pf32  "  "*2^23  "n23  (B157) 

To  write  the  moment  equation,  let 


where 


Then, 


•• 


N32* 


N. 


P3 


N, 


G2 


(B158) 


(B159) 


F23&b2 


W82F23a2  -  F12rb2  +  W8iFI2(d1  “  al>  +  W>2('*2  +  V 

h*3l* 


(B160) 


The  conservative  bearing  forces  F  -  and  F  .  are  obtained  from  the  simulta¬ 
neous  solution  of  the  following  component  equations  of  the  force  expression  B153 
for  Fx2  and  Fy2: 


-F23sin(B2  +  82)  -  F23us2cos(B2  +  0^  -  F^sin  ((^  -  0  )  (B161) 

-♦-us iF i2c°s C 3 1  -  0j)  +  Fx2  +  uFy2  +  T2cosy2  -  0 

F23co®( e2  +  V  "  U82F238inC62  +  V  +  Fi2cos(0l  ~  V  (B162) 
+p8lFl28in<6l  “  V  +  Fy 2  -  uFy2  +  T2siny2  -  0 

'  This  leads  to: 

^  Fx2  "  F?3A36  +  F12A37  +  T2A38  (B163) 
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Fy2  "  F23A39  +  F12A40  +  T2A41  (B164) 

^A^f  *irfn*  *r®  Positive  in  order  to  furnish  positive  friction  moments  in 

equation  B160*)  In  the  above: 

2 

,  C I  -  M  a2)sin(80  +  62)  +  u(l  +  s0)cos(B,  +  8„) 

A36  “  I  “ - 2 - - - -  I  03165) 

I  +  u 

(  (1  +  w^sin^  -  0,)  +  u(l  -  s,  >00^(8,  -  0,) 

A37  "  1  “ - 2 - - - -  I  (B166) 

1  +  U 


COSY 2  -  ysiny  2 

A38  "  *  “  j  +  '2  I  (B167 ) 

A  ,  +  82)8ln(fil2  +  92^  “  U  “  U282)cos(B0  +  eo) 

39  "  I  - - 2 - -  I  (B168) 

1  +  u 

a  i  ”  8P8in(si  ~  8j)  ~  (1  +  u2s , )cos(&,  -  of) 

A40  “  I  —"2  - - - —  I  (B169) 

1  +  M 


UCOsy2  +  siny2 


1  +  p 


2 


(B170) 


Substitution  of  equations  B163  and  B164  into  the 
subsequent  solution  for  the  contact  force  F23  gives 


23 


F12A42  "  T2A43  +  I2N32<^ 
^ - 


moment  equation  B160  and 


(B 171) 


where 


*42  "  rb2  ‘  u(’l(dl  *  al>  +  l>2<A37  *  a40)I  (B172) 
*43  ■  w2  (*38  +  *41>  (B173) 
*44  "  “b2  *  ul*2*2  *  °2(*36  +  *39»  (B174) 
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Dynamic*  of  Combin«d  System  in  Coupled  Motion  (applicable  to  both  conflguratione) 


To  obtain  a  single  differential  equation  for  the  totel  system  In  coupled 
motlou  in  terms  of  the  escape  wheel  engle  an  appropriate  expression  for  the 
contact  force  F,*,  which  aloo  contains  the  contribution  of  the  rotor,  must  be 
substituted  into  either  equation  B100  or  equation  B119  for  entrance  or  exit- 
coupled  motion,  respectively. 

Thus,  first  substitute  equation  B144  for  F^2  Into  equation  B171,  which  is 
the  above  expression  for  Fgg: 


23 


A42A35 


“A42A32  „  .  A42A33  ?  A42A34  „2  r2 

“ST '  2  43  '  “  ~  31  * 


*3,  sln(*lc  -  N3X<>T>  +  «2*32  '  7^  I1RN31)*  1 


(B175) 


Now,  equation  B175  is  substituted  into  equation  B100  or  equation  B119  and  the 
final  differential  equation  of  coupled  motion  results: 


a455  +  +  A47* 

"  A48  “  A49flin<^lc  +  N31 V  +  A508in(Yi  "  *  “  V 


where 


B-5 


(B176) 


A45  "  X3A18  +  XPRA17U  "  A^4  (I2N32  “  A3J  I1RN31) 


(B177 ) 


_  A, _A. „A_. A, „  _ 

. 2  .  15  18  34  42  „2 

A , ,  *  InuA.-jV  +  A.,A01U  +  1  . . . —  N,, 

46  PR  17  17  21  A^A^  31 


(B178) 


,2  A  A  A  A 

A  u  2  r^~r  A  A  U  +  —  ~  ^3  Jiik  M 

A47  1  TSTf  A17A20U  +  A31A44  31 


(B179) 


The  coefficients  A4^  to  A^q  are  only  valid  for  entrance  contact;  i.e.,  for 
equation  B100.  If  exit-coupled  motion  is  in  operation,  equation  B 1 1 9  is 
applicable  and  k^-j  becomes  AA^,  and  Ajg  becomes  AAjg. 
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(B180) 


*  ~A15A18A32A42  A13A18A43  ,  A  .  T 

A48  A44  T2  "  A16A18T3 


u2ai7a19 


,  A15A18A35A42 

49  A31A44 


(  B 1 8 1 ) 


m  r  fl.A.-w 
p  cp  4  17 


(B182) 


Contact  Force  Equations  for  Coupled  Motion  (applicable  to  both  configurations) 


The  contact  force  F2o  is  given  by  equation  B175: 


23 

-  _L_  [  ~A32A42 


F23  A44 

A35A42 


'31 


A33A42  •  A42A34  2  *2 

-  T2A43  -  “TT  N31*  "  "XT  N31* 


'31 


'31 


'42 


A31  8in(*lc  +  N3lV  +  (I2N32  A31  *1R"31 


(B183) 


The  contact  force  F^  is  found  with  the  help  of  equation  B171: 

•• 

F23A44  +  T2A43  “  12N32<*> 


F  m 
*12 


'42 


(B184) 


The  contact  force  Pn  is  given  either  by  equation  B96  or  equation  B117: 


Pn  " 


-V  +  F23A15  ~  Vl6 
A17 


(B185) 


Note  that  for  exit-coupled  motion;  i.e.,  equation  B117,  the  term  A^7  becomes 
AA^y.  If  this  force  is  desired  in  terms  of  the  pallet  variable  \j>,  equation  B76 
or  equation  B107  must  be  used: 


V  *  A2i*2  +  2  TCT 


+  uTa,  „  -  m_r 


'20’ 


‘19 


JL£L 


4“  "  *  -  *c) 


'18 


Again,  note  that  for  exit  motion  Ajg  becomes  AAjg. 


(B186) 


■*” VTIW.  T  ~  nw'  *” 


Differential  Equations  for  Free  Motion  (applicable  to  both  configurations) 


Pallet  Free  Motion 


By  letting  Pn  -  0  in  equation  B186,  the  single  free  motion  differential 
equation  for  the  pallet,  which  is  now  Independent  of  exit  or  entrance  contact,  is 
obtained* 


a5]  J  +  A21$2  +  A52l  -  -A53  +  A54sin(Y^ 

where 

a51  "  XPR 

A52  "  2  T«r  A20 
A53  "  ^19 
A54  *  “Vcp41* 


(B187) 

(B188) 

(B189) 

(B190) 

(B191) 


Escape  Wheel  -  Gear  Train  -  Rotor  Free  Motion 


By  setting  Pn  ■  0  in  equation  B96  or  equation  B117,  the  common  differen¬ 
tial  equation  of  the  system  without  the  pallet  is  obtained  (again  Independent  of 
entrance  or  exit  contact).  This  leads  to: 


V  ’  A15F23  "  T3A16 


(B192) 


Equation  B175  is  then  substituted  into  the  above  for  F23.  The  resulting 
differential  equation  is  given  by: 


A55*  +  A56*  +  A57*  ‘  A58  '  A598ln(*lc  +  “31V 


(B193) 


where 


A15  A15A42 

Ac c  -  I-  -  t**  N„I-  +  T-\  N..  1 

55  3  A44  32  2  Aj^A44  31  1R 


(B 1 94 ) 


(B195) 


,  A15A34A42  n2 

56  A31A44  31 


‘57 


A15A33A42 

A31A44 


N 


31 


A58  "  "A16T3  ' 


A  A 
15  43 

A44 


t  - 
2 


A15A32A42 

A31^44 


‘55 


A15A35A42 

— T — T - 

A31A44 


(B196) 


(B197) 


(B198) 


Contact  Force  Equations  for  Free  Motion  (the  subscript  F  stands  for  free  motion) 
(applicable  to  both  configurations) 


Equation  B192  may  be  solved  for  the  free  motion  contact  force  Fp23i  once  <J> 
is  known : 


FF23 


V  +  T3A16 


(B199) 


Equation  B171, 
and  pinion  set  no. 

FF12: 


FF12 


which  was  derived  from  the  force  and  moment  equations  of  gear 
2,  may  be  modified  to  obtain  the  free  motion  contact  force 


FF23A44  +  T2A43  “  12*32* 


'42 


(B200) 


It  must  be  understood  that  in  both  equations  B199  and  B200  the  angular  accelera- 

•• 

tlon  of  free  motion  <j>  as  obtained  from  the  solution  of  equation  B193  is  used. 


Impact 


While  the  following  will  show  that  the  general  form  of  the  Impact  equation 
is  the  same  for  entrance  and  exit  contact,  it  must  be  kept  in  mind  that  the 
expressions  for  the  moment  arms  Aj,  Bj,  C!  and  D!  depend  on  whether  the  angle  a 
or  the  angle  «EX  is  Involved.  (See  equations  D4,1  D5,  D8,  and  09  in  appendix  0.; 
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Entrance  Impact 


The  general  condition  of  entrance  Impact  le  Illustrated  In  figure  B- 
11.  As  In  reference  B-2,  It  Is  assumed  that  the  only  acting  Impulse  on  the 
pallet  and  the  escape  wheel  Is  represented  by  the  mutual  normal  Impulse  P 
between  these  mechanism  components.  Experience  has  shown  that  the  frictional 
Impulse  uP  may  be  disregarded.  Under  these  circumstances ,  the  angular  impulse 
on  the  pallet  becomes: 

Jn  -  Djn  x  P  n  -  P  Djk  (B201) 

p  It  n  n  n  1 

The  angular  Impulse  on  the  escape  wheel  Is  given  by: 

J.  -  AX  x  (-P  )n  -  -P  A'ic  (B202) 

3  It  n  u  n  1 

If  the  above  equations  are  made  part  of  the  usual  angular  momentum  relationships, 
the  components  can  be  expressed  as 

Ip($f  -  ♦i)  -  PnD]  (B203) 

and 

I3(Jf  -  -  -PnA{  (B204) 

where  I  and  1^  represent  the  polar  mass  moments  of  inertia  of  the  pallet  and  the 
escape  Vheel,  respectively.  The  subscript  1  indicates  the  angular  velocities 
before  impact  while  f  indicates  the  same  quantities  after  impact.  When  Pn  is 
simultaneously  eliminated  from  equations  B203  and  B204,  the  following  expression 
results : 


va;*£  +  Wf 


w +  r3D;+i 


(B205) 


Exit  Impact 


The  conditions  of  exit  impact  are  shown  in  figure  B-12.  With  the  same 
assumptions  as  in  part  a  above,  the  expression  for  the  angular  impulse  on  the 
pallet  becomes: 

j  -  Djn  x  (-P  )n  -  -P  D'k  (B206) 

p  It  n  n  n  1 
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usstgc wvmtsTP »gai^g^r ww.'rn  wj. '«;•"%■■•-  i.n*n vw iw* 


The  angular  Impulse  on  the  escape  wheel  now  becomes: 


J3  -  *Jne  *  (Pn)nn  -  Pn*Jk 


(B207) 


Again,  these  expressions  are  used  In  the  angular  momentum  formulations  of  both 
components,  as  follows: 


y** 


-p  d; 

n  1 


(B208) 


*3(*f 


(B209) 


When  Pn  is  now  simultaneously  eliminated  from  equations  B208  and  B209,  equation 
B205  results  again.  Thus,  Identical  formulations  may  be  used  for  entrance  and 
exit  contact. 


Impact  Equations  for  Escapement 


.  t  In  order  to  obtain  expressions  for  the  post  impact  angular  velocities 

4>f  and  ipf,  the  restitution  equation 


V  -  V 

TNf  SNf 

't  V  -  V 

vTNi  vSNi 


(B210) 


must  be  solved  simultaneously  with  equation  B205*  Vm™  and  V-™  are  defined  by 
equations  DI3  and  D15  in  appendix  D.  The  velocities  VTNf  and  VgNf  represent  the 
normal  contact  point  velocities  after  impact.  After  appropriate  substitutions, 
the  simultaneous  solution  of  these  expressions  furnishes: 


♦l(I3°;2  -  *rIpA{2)  +  ♦lI.AI(1  +  eg) 
VI2  ♦  Vi2 


(B211) 


*fAi  -  er(^Di  -  ^Aj) 


(B212) 
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Changes  in  Impact:  Equation  due  to  Presence  of  the  Rotor  and  Gear  and  Pinion 
No.  2 


The  presence  of  the  rotor  and  gear  and  pinion  no.  2  is  accounted  for  by 
referring  their  moments  of  inertia  to  the  escape  wheel  shaft.  Thus,  the  total 
escape  wheel  moment  of  inertia  is  given  by: 

IST0T  ■  +  vL  +  V31  (B213> 
Equations  B122  and  B159  contain  the  above  transmission  ratios. 
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Nomenclature 


The  following  derivations  for  the  kinematics  of  coupled  motion  are 
applicable  to  both  entrance  and  exit  motion.  It  must  be  remembered  that  for 
entrance  motion,  according  to  figure  C-l: 


a  *  a 

en 

(Cl) 

and 

S  “  Sen 

(C2) 

while 

for  exit  motion: 

a  **  a 

ex 

(C3) 

and 

8  "  «ex 

(C4) 

Further,  from  figure  C-l 

<|>  -  Escape  wheel  angle,  measured  with  respect  to  positive  x'  -  axis. 

-  Verge  angle,  measured  with  respect  to  positive  x'  -  axis.  For 

entrance  condition,  defined  by  line  0  V.  For  exit  condition, 
defined  by  line  0  W  (measured  ccw  with  respect  to  the  positive  x'  - 
axis).  p 

a  *  Distance  between  pivot  points  0p  and  08 

b  -  Escape  wheel  radius 

c  *  For  entrance  condition,  length  0pV.  For  exit  condition,  length  0  W 
(usually  identical).  p 


Unit  Vectors 


Tangential  ard  normal  unit  vectors,  which  are  associated  both  with  the 
entrance  and  the  exit  working  surfaces,  are  shown  in  figure  C-l.  When  expressed 
without  distinction  between  these  surfaces,  they  become: 


n£  »  cos(ip  +  «)Z*  +  ein(ip  +  a)J ’ 


and 


n  -  -sin(ip  +  a)!'  +  cos(ip  +  a)J ' 

D 


(C5) 

(C6) 


95 


(C6) 


n^  ■  — ain(  ^  +  a)i*  +  cos(tp  +  a)J* 

Further,  the  unit  vectors  In  directions  08S  end  0  V  (or  0  W),  respectively,  are 
given  by  P  P 


and 


»  cos  <p  i*  +  sin  <J>  j* 

n  ■  cos  ip  i'  +  sin  ip  J* 
c 


(C7) 

(C8) 


Loop  Equation  and  Input-*Output  Relationship 


The  loop  equation  for  coupled  motion  is  given  below.  Contact  point  S  Is  at 
the  tip  of  the  escape  wheel  tooth.  The  coincident  point  on  the  pallet  is  point  T 
(fig  C-l). 


OV+VT+SO  +00  -0  (C9) 

p  S  8  P 

Substitution  of  the  appropriate  dimensions  and  unit  vectors  into  the  above 
furnishes: 


cn  +  gn.  -  bn,  -  ai*  •  0  (CIO) 

C  t  D 


This  becomes  with  the  help  of  equations  C5  to  C8: 


c(cosipI'  +  sinipj')  +  g(cos(ip  +  o)i'  +  sin(ip  +  a)J')  (Cl  1 ) 

^(cosipi'  +  sinipj*)  -  ai'  ■  0 

Appropriate  ordering  of  equation  Cll  leads  to  component  expressions  which  may  be 
solved  for  g  and  ip: 


and 


bcosip  +  a  -  ccosip  -  geos  (ip  +  a)  -  0 


bsinip  -  csinip  -  gsin(ip  +  a)  ■  0 


(C12) 

(C13) 
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Both  of  those  component  equations  may  be  aolved  for  g.  Thue, 


.  bcoa$  +  a  -  ccoa\ 
®  "  cosTT  +  a) 


(C14) 


bsinA  -  csinii 


(C15) 


g  may  be  eliminated  by  setting  equation  C14  and  C15  equal  to  each  other.  This 
leads  to 

(bcos$  +  a  -  ccostJ))sin(tJ>  +  ot)  ■  (bsin<|>  -  csini|>)co8(ij;  +  a)  (C16) 

After  expansion  and  rearrangement,  a  four-bar  linkage  type  expression  is  obtained 
which  may  be  solved  for  the  angle  according  to  the  method  shown  by  R. 
Hartenberg  and  J.  Denavit.* 


AsiniJ/  +  Bcosip  ■  C 


(C17) 


where 


A  -  a  cosat  +  b  cos(<|>  -  a) 
B  ■  a  sincx  -  b  sln(<j>  -  a) 


C  ■  c  sina 


The  solution  to  equation  Cl 7  is  of  the  form: 


<J>  ■  2tan 


-1  A  1  /A2  +  B2  -  C2 


(C18) 


Once  the  correct  <p  has  been  obtained  for  a  given  a  or  a ,  the  associated  gen 
or  gex  may  be  found  with  the  help  of  equation  C14  o®nequatfon  C15.  For  the  given 
directions  of  the  associated  unit  vectors  gfin  will  be  a  positive  quantity  while 


gex  will  be  a  negative  one. 


*  Kinematic  Synthesis  of  Linkages,  McGraw-Hill  Publishing  Co.,  New  York,  1964. 
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Angular  Velocity  of  Pallet 


Implicit  differentiation  of  equation  C17  leads  to  the  following  expression 
for  the  angular  velocity  of  the  pallet  when  driven  in  coupled  motion  by  the 
escape  wheel: 


where 


•  •  P 

♦  ■  ♦  ^ 


■  Angular  velocity  of  the  escape  wheel 
P  ■  bcos(4>  -  a  -  y) 

Q  ■  acos(i|>  +  a)  +  bcos($  -  a  - 


(C19) 


(C20) 

(C21) 


Relative  Velocity  of  Contact  Point  S  on  the  Escape  Wheel  with  Respect  to  Contact 
Point  T  on  the  Pallet 


With  point  S  the  tip  of  the  escape  wheel  tooth  and  point  T  the  coincident 
contact  point  of  the  pallet,  the  relative  velocity  Vg/T  1°  the  direction  of  the 
unit  vector  nt  is  found  from  (fig.  C-l): 

^S/T  "  [VS/08  *  nt  "  ^T/0p  *  nt]  nt  (C22) 

The  velocity  Vg/Q  of  point  S  with  respect  to  the  escape  wheel  pivot  is 
obtained  from  8 


vs/o  ■  +k  *  bnb 


(C23) 


while  VT/0  ,  the  velocity  of  point  T  with  respect  to  the  pallet  pivot,  is  given 
by:  P 


VT/0  '  ,E  *  <C"c  +  S"t) 
P 


(C24) 


Substitution  of  equations  C23  and  C24,  together  with  equations  C5,  C7,  C8,  into 
C22  leads,  after  the  appropriate  vector  manipulations,  to: 


VS/T  "  4»bB±nC  4»  -  ♦  -  a)  +  ^csinajnt 


(C25) 


WKVnVWUUV!!-'  ;TTI'.’T^^Vf'T1\TTlT 


y~v*ff  t  w£7.  rr^7 — vj’rKj’-'«?;7 -;w--:’«^' 


i 


I 


y 


s 


ft-. 


It  can  be  shown  that  the  above  expression  also  represents  the  vector  g  - 
gn  ,  as  obtained  from  the  rate  of  change  of  the  quantity  g»  which  was  defined  by 
aquations  C14  and  CIS* 

Angular  Acceleration  of  the  Pallet 

Differentiation  of  equation  C19  with  respect  to  time  results  in  the 
following  expression  for  the  pallet  angular  acceleration: 


••  •• 


1|/  -  u*  +  v$2 


where,  with  equations  C20  and  C21 


(C26) 


and 


U  -7T 


V  ■  ~-r  [aP2sin(<i)  +  a) 


-  b(P  -  Q)  sin(<(>  -  a  -  0')] 


(C27) 


(C28) 


u 

$ 

K 


tv 

f." 


KW 


l*: 
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Moment  Arms 


The  following  derivations  for  the  moment  arm  vectors,  which  find  use  in  both 
coupled  motion  and  impact  calculations,  are  valid  for  both  entrance  and  exit 
conditions  (figs.  D-l  and  D-2)  as  long  as  the  appropriate  angle  n  or  n  is 
used.  The  variables  <j>  and  <{/  must  be  available  from  the  free  motion  computations. 

The  following  loop  for  the  determination  of  vectors  and  Bj  is  shown  by 
figure  D-l ; 


A'n  +  B*n  -  bn.  -  0  (Dl) 
it  in  o 

Substitution  of  the  unit  vectors  given  by  equations  C5  through  C7  and  subsequent 
separation  of  component  expressions  leads  to  the  following  two  equations  which 
may  be  solved  for  AJ  and  Bj: 


AJcos(iJ)  +  a)  -  B|sin(>|)  +  a)  ■*  b  cos  <f> 


and 


Ajsin(i^  +  a)  +  Bjcos(i{>  +  a)  -  b  sin 
Simultaneous  solution  of  the  above  expressions  for  A|  and  Bj  leads  to 
AJ  -  bco8((|>  -  \p  -  a)nt 


and 


BJ  -  bsin(<j>  -  t|>  -  «)nn 

Figure  D-l  also  serves  for  the  determination  of  the  vectors  Cj  and  D! . 
equation 


(D2) 


(D3) 


(D4) 

(D5) 


The  loop 


Dlnt  +  Clnn  ~  cnc  -  gnt  -  0 

furnishes,  with  the  help  of  equations  C5,  C6,  and  C8  for  the  unit  vectors,  the 
following  set  of  component  expressions: 


(DJ  -  g)cos(ij>  +  a)  -  CJsin(if>  +  a)  ■  c  cos 


(Db) 


and 
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(D|  -  g)sln(ij/  +  a)  +  C|cos(\|»  +  a)  -  c  sin 

(D7) 

Simultaneous 

solution  for  C{  and  D[  gives: 

C*  »  -csinon 

1  n 

(D8) 

and 

Dj  ■  (ccosa  +  g)nfc 

(D9) 

Again,  It  must  be  remembered  to  substitute  the  appropriate  entrance  or  exit 
values  for  a  and  g. 


Normal  Velocities  of  Contact  Points  Before  and  After  Impact 


The  contact  point  S  on  the  escape  wheel  and  the  contact  point  T  on  the 
pallet  are  shown  In  figure  I>-3.  At  the  Instant  before  Impact,  the  respective 
escape  wheel  and  pallet  angular  velocities  are  given  by: 


-  (j^k’  (DIO) 

and 

-  ^k'  (Dll) 

The  associated  velocity  component  Vgjj^  of  the  escape  wheel  contact  point  S  Is 
obtained  from: 


^SNi  "  *1  X  (D12) 

This  velocity  component  Is  normal  to  the  pallet  face.  This  becomes,  with  the 
help  of  equations  DIO  and  D4: 


VgNi  -  b^cosO  -  i|»  -  a)nn  (D13) 

In  a  similar  manner  the  velocity  component  of  the  pallet  contact  point  T 

becomes : 


(D14) 
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This  velocity  component  is  also  normal  to  the  pallet  face.  Substitution  of 
equations  D9  and  Oil  into  the  above  furnishes: 


Vtn1  “  ^(ccosn  +  g)nn  (015) 

When  equations  D13  and  D15  are  adapted  to  the  ^post-impact  normal  velocities 
^SNf  an<*  ^TNf*  t*ie  respective  angular  velocities  ^  and  must  be  substituted 
(impact  section  of  appendix  B). 


The  contact  sensing  quantities  f  and  g'  for  entrance  and  exit  conditions  are 
defined  by  figures  E-l  and  E-2,  respectively.  It  is  assumed  that  the  angles 
$  and  are  known  at  every  instant  of  free  motion.  To  determine  the  entrance 
side  sensing  quantities,  the  following  loop  equation  is  written: 

+  "ST  +  W  +  W  +  0  0  -  0  (El) 

s  p  p  s 

This  becomes  with  the  help  of  the  appropriate  unit  vectors 

bn.  +  fn  -  g'n,.  -  cn  +  ai'  -  0  (E2) 

b  n  t  c 

When  a  similar  expression  is  written  for  the  exit  side  (fig.  E-2),  the  following 
relationship  is  obtained: 

bn.  -  fn  +  g'n,.  -  cn  +  ai'  ■  0  (E3) 

b  n  °  t  c 

The  quantities  f  and  g'  have  opposite  signs  in  equations  E2  and  E3.  If  equation 
E2  is  solved  for  f  and  g'  with  the  help  of  the  unit  vectors  of  equations  C5  to 
C8,  then  for  the  entrance  side: 


f  -  asin(ij>  +  a)  -  bsin($  -  ip  -  a)  -  csina  (E4) 

and 

g'  -  acos(i^  +  a)  +  bcos(<|>  -  \p  -  a)  -  ccosa  (E5) 


The  above  expressions  will  yield  positive  numerical  quantities.  When  the*  same 
expressions  are  used  for  the  exit  side  (with  the  appropriate  angle  a  ),  the 
resulting  numerical  quantities  will  be  negative  because  of  the  aforementioned 
sign  differences  between  equations  E2  and  E3. 


Ill 


rrTV'vp, 


Angles  A,  a  ,  end  a 
en  ex 

The  noraencleCure  of  the  lndlceted  straight  sided  verge  is  defined  in  figure 
F-i.  The  angle  A  is  obtained  as  follows: 


A  -  2tan 


-1  h 


With  the  above,  and  angle  n,  the  angles  otfln  and  become,  respectively: 


a  ■  TT  -  -r- 

en  2 


“e*  "  7"  n 


The  angle  n  is  a  basic  design  dimension. 


Lengths  of  Working  Surfaces 


Figure  F-2  is  used  to  find  the  distance  pen;  i.e.,  the  length  of  the 
straight  portion  of  the  entrance  working  surface,  as  well  as  the  corresponding 
length  p_v  of  the  exit  working  surface  of  the  verge. 

'“X 

The  length  p  -  LV  is  determined  from  the  difference  in  lengths 


where 


PeB  ■  y(L')  -  f 


y(L')  -  the  y-coordinate  of  the  point 
yl  ■  r  ,  and  the  straight  line 


sint  of  intersection  of  the  circle  x4  + 
line  x  »  L’ 


Then: 


L"2  +  y2  -  r2 


'  2  2 
r4  -  V4 


Equation  F4  becomes: 
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UI 


EXIT  WORKING  SURFACE 


".n  - -  L'2  -  T  <"> 

The  length  pQX  ■  MW  It  obtained  with  the  help  of 

Pex  -  W»  -  t  (F6) 

where  point  M  represents  the  blend  point  between  the  fillet  of  radius  q  and  the 
straight  line  of  the  exit  working  surface.  The  length  W'  is  obtained  from 

W  -  (F7) 

cosn 

while  t  is  given  by: 


IT  , 

j  +  n 

tan  ( - j— ) 

where  |  +  n  is  the  angle  spanned  by  the  fillet, 
becomes 


With  the  above,  equation  F6 


L'  -  m 


f+n 

tan  ( — a — ) 
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